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Nanomechanical oscillators have been employed as transducers to measure force, mass 
and charge with high sensitivity. They are also used in opto- or electromechanical 
experiments with the goal of quantum control and phenomena of mechanical systems. 
Here, we report the realization and operation of a hybrid monolithically integrated 
transducer system consisting of a high-Q nanomechanical oscillator with modes in the 
MHz regime coupled to the near-field of a high-Q optical whispering-gallery-mode 
microresonator. The transducer system enables a sensitive resolution of the nanome- 
chanical beam's thermal motion with a signal-to-noise of five orders of magnitude and 
has a force sensitivity of 74aNHz^^^^ at room temperature. We show, both theoreti- 
cally and experimentally, that the sensitivity of continuous incoherent force detection 
improves only with the fourth root of the averaging time. Using dissipative feedback 
based on radiation pressure enabled control, we explicitly demonstrate by detecting a 
weak incoherent force that this constraint can be significantly relaxed. We achieve a 
more than 30-fold reduction in averaging time with our hybrid transducer and are able 
to detect an incoherent force having a force spectral density as small as ISaNHz^^''^ 
within 35 s of averaging. This corresponds to a signal which is 25 times smaller than 
the thermal noise and would otherwise remain out of reach. The reported monolithic 
platform is an enabling step towards hybrid nanomechanical transducers relying on the 
light-mechanics interface. 



Nanomechanical oscillators [T] serve as ultrasensitive de- 
tectors of force [2j, mass and charge [1]. Recently, in- 
creasing efforts have been devoted to sensitively detect the 
nanomechanical motion of these oscillators [SHS] with recent 
systems exhibiting a sensitivity below that at the standard 
quantum limit (SQL) O [lOj. For sensitive force detection 
the requirements on displacement sensitivity are less strin- 
gent, though, because of the thermal limit. Recent work 
has demonstrated that trapped ions have also the potential 
to be employed as sensitive transducers for ultrasmall forces 
[111 I12| . with the force sensitivity reaching levels of only 
SyNHz"^''^ [12]. While being a promising approach to de- 
tect specific small forces, it presently suffers from challenging 
technology required for ion trap experiments, stable interac- 
tion times being only in the millisecond range [TT] and a lack 
of field-deployable sensors. In contrast, cantilever-based sens- 
ing is a well-established technique that allowed remarkable 
achievements such as the detection of single spins [O] and 
the reconstruction of the structure of a virus [14], and it is 
increasingly used in the emerging field of biosensing |15| . A 
particularly promising approach is to parametrically couple a 
nanomechanical oscillator to a microwave |9|[T6] or optical [17j 
microcavity, thus enhancing the displacement sensitivity and 
allowing to efficiently resolve the motion of cantilevers with 
dimensions below the diffraction limit. Previous realizations 
suffered, however, from the necessity of cryogenic operation 
[21 [TB] , from the fact that the optical resonator and mechani- 
cal oscillator were realized on two different chips and had to be 
brought to close proximity by external positioners [17] mak- 
ing the system too unstable for usage as a force transducer 
or from low mechanical Qm, significantly limiting the force 
resolution |18l I19j . Here, we present an integrated hybrid 
transducer system consisting of a high-stress Si3N4 nanome- 
chanical beam whose high-QM out-of-plane mode is coupled 
to the near-field of a whispering gallery mode of a silica disk 



resonator. As such, we combine the arguably best material 
for high-QM nanomechanics [2D] (stoichiometric Si3N4) with 
an advantageous material for linear optics [21H24| (Si02) in 
an integrated on-chip device, which allows the detection of 
ultra-small forces at the aN level. The advantages of silica for 
optical resonators result from higher thresholds for thermal 
nonlinearities when compared to silicon f2F or silicon nitride 
[26| as well as the absence of two photon absorption JJ, i thus 
allowing significantly higher circulating power and preventing 
cross talk between optical modes induced by the nonlinear- 
ity. Moreover, silica microresonators can be coupled with high 
efficiency to optical fiber as required for highly efficient mea- 
surements of position or force. 

The versatility of the presented hybrid transducer is 
enhanced by implementing a feedback scheme that intro- 
duces an effective susceptibility for the mechanical resonator 
[28H30] . This enables tuning of the mechanical resonance 
frequency and linewidth. This feature is particularly in- 
teresting in the context of ultra-sensitive force detection, 
which relies on the use of high-Qivi mechanical oscillators 
whose detection bandwidth is limited to a very narrow 
range around their resonance frequency [281 I31| . Here, we 
report the implementation of an integrated feedback scheme 
enabling full and independent control of the effective me- 
chanical properties of the transducer. Though being limited 
mostly by thermal noise, nanomechanical transducers can 
be used to detect arbitrarily weak coherent signals, with the 
sensitivity scaling with the square root of the averaging time 
T [13] . The situation is more delicate for detecting incoherent 
signals, such as quantum backaction T^, atomic mass shot 
noise [3], or a gravitational wave background [33[ . Here the 
absence of phase coherence for such signals prevents using 
amplitude averaging to retrieve their signature, therefore 
energy averaging must be employed [13]. However, this is 
made at the expense of the measurement resolution, which 
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increases only very slowly with t^^^ |13) . Here, we propose 
and experimentally demonstrate with our hybrid transducer 
that dissipative feedback allows a significant decrease of 
the averaging time required to reach a given accuracy 
in the context of a low signal-to-noise ratio (SNR) force 
measurement. 

Near-field coupled hybrid transducer. Our integrated 
system enables hybrid integration of the approach of Anets- 
berger et al. [T7] and is fabricated using a simple procedure 
easy to be implemented for mass production (fabrication de- 
tails are given in Methods). Images of the structure are given 
in Figures 1(a) and 1(b). It consists of a doubly-clamped, 
high-stress Si3N4 nanomechanical beam 120j (blue) which is 
realized above the wedge of a disk made of Si02 (red) . Both, 
the pads to which the oscillator is clamped and the silica disk 
rest on silicon pedestals (green). We choose this geometry to 
optimize the hybrid transducer's performance as a force sen- 
sor, which is directly related to the mechanical quality factor 
Qm and mass of the cantilever. It was shown [50] that the fun- 
damental out-of-plane mode of a Si3N4 nano-beam exhibits 
distinctly high Qm factors while having a small mass, thus 
making it an optimal candidate for a force transducer. In or- 
der to efficiently couple the out-of-plane motion of the beam to 
the evanescent field of a WGM disk cavity, the beam needs to 
be positioned above the wedge of the disk resonator, where the 
distribution of the near-field electromagnetic energy is maxi- 
mized (cf. Fig. 1(d)). The beam of the hybrid system used in 
this work has dimensions of 90 /im x 700 nm x 100 nm and the 
disk has a diameter of 76 /im, with the distance between the 
beam and the disk being approximately 250 nm. Like micro- 
toroidal cavities [23j silica microdisk cavities support high-Q 
optical whispering gallery modes (WGMs) [34] . In our system 
we measured different families of WGMs around a wavelength 
of 1550 nm by coupling laser light coming from an external 
cavity diode laser to the cavity via the fiber-taper technique 
|35| (Fig. 1 (c)). The linewidth of such modes is around 
Ji/ « 3 GHz corresponding to a value of Q = v/Si/o ~ 65, 000 
with vo being the resonance frequency of the mode. This 
measured linewidth is considerably higher when compared to 
previous results |24l I34| , which is explained by enhanced scat- 
tering losses due to the defects present on the disc arising from 
shadowing effects during fabrication. While the increased 
optical linewidth limits the attained displacement sensitiv- 
ity of the mechanical motion to the level of several tens of 
fmHz"'^''^, the force sensitivity is not affected, since the ther- 
mally driven motion of the mechanical transducer is well re- 
solved as shown below. In fact, an improved linewidth would 
even be detrimental for the system's ability to function as a 
force transducer due to effects of dynamical backaction [36], 
which would affect the stability of the mechanical transducer. 

To measure the transduction of the mechanical motion of 
the oscillator, we keep the hybrid system under a pressure 
of below 10"'' mbar and detune the laser to the blue side 
of the fringe of a critically coupled optical mode [35]. We 
are able to resolve thermally driven mechanical motion 
of the nanomechanical oscillator, which is imprinted as a 
modulation on the transmitted optical intensity. Launching 
a laser power of 400 into the cavity, we observe the 
fundamental out-of-plane mode of the beam at a resonance 
frequency of JIm/Stt = 2.88 MHz with a signal-to-noise ratio 
of 50 dB. By applying a known phase modulation ^37) we 
can determine the optomechanical coupling of the system 
to be G/2n — dy/dx = 2.9 MHz/nm (x is the displacement 



of the mechanical oscillator) and calibrate the spectrum 
in frequency noise Sw [^] as well as displacement density 
S^rc [^] = Sr,^ [r2] /G^ units. The eflFective mass [3^ rricfi of 
the fundamental out-of-plane mechanical mode coupled to 
a WGM optical mode equals half the physical mass m |17) . 
thus giving mefi = 9 x 10"'^ kg for our system. By fitting 
the data we obtain a very high mechanical quality factor 
of Qm = 4.3 X 10^. To confirm this value, we perform a 
ringdown measurement of the oscillator by actuating the 
beam at JIm via the radiation pressure force and interrupting 
the actuation while recording the response (scheme is shown 
in Fig. 2(a)). We obtain a value of 4.8 x 10^ (Fig. 2(b)) 
corroborating the spectrally obtained result. To exclude 
dynamical backaction effects via radiation pressure in our 
measurements, we varied laser power and detuning around 
their working points and did not observe any change in Q.m 
or Qm, ensuring that this type of interaction is indeed weak. 

Actuation via the radiation pressure force. To 

demonstrate the actuation of the mechanical oscillator via 
the radiation pressure force, we perform a force response 
measurement using a pump-probe scheme. The setup is 
shown in Fig. 2(a) and details of the experiment are given 
in Methods. Fig. 2(c) shows the response of the mechanical 
oscillator to different levels of actuation. The response 
increases with stronger actuation and shows a nonlinear 
behavior above a given threshold that can be described by 
the Dufling model [39| . 

Feedback tuning and cooling of the transducer. We im- 
plement a feedback scheme i29^ to effectively address the me- 
chanical properties of our nanomechanical oscillator, which, 
as shown later, is important for an improved force detection. 
Feedback consists in applying a force Ffb to the oscillator 
proportional to its displacement x. In Fourier space this force 
can be written as: 

Ffb[n] = -ye** X meftTMfiMa;[n], (1) 

where g and $ denote the amplitude and the phase of the gain, 
respectively. Thus, in presence of feedback, the response of 
the oscillator to an external force _Fcxt reads: 

x[n] (-Fcxt [!^]-(-Ftb[n]) 

^X^«[^]F.A^l (2) 

where Xcff denotes the (Lorentzian) effective susceptibility 
[40] . characterized by effective mechanical resonance fre- 
quency r^cff and effective damping given by: 

n,\ = r!^fl + -^cos<&') (3) 
\ Qm / 

Fcff =rM(l+.gsin$). (4) 

The real part of the gain g cos $ determines the conservative 
contribution of the feedback force modifying the effective me- 
chanical resonance frequency of the oscillator as shown by 
Equation [S] Equivalently, the imaginary part of the gain 
(;sin<I> sets the dissipative component of the feedback force 
changing the effective mechanical damping rate of the oscilla- 
tor (Eq. |4]|. Consequently, independent control of amplitude 
g and phase $ of the feedback enables tuning the oscillator's 
mechanical response. 
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Historically feedback implementation has been based on 
the use of resonant electrical circuits designed to create a feed- 
back signal at the mechanical frequency of interest [3UII41lB^ . 
This method suffers from the difficulty to tune these circuits, 
so that the actual feedback implementation is limited to one 
specific mechanical mode. Here, we use a method based on a 
demodulation-modulation technique [4:3 , which enables accu- 
rate control of any mechanical mode up to arbitrary high fre- 
quencies. The scheme employed is shown in Figure 3(b). The 
method relies on reading out the time evolution of the quadra- 
tures Xi and X2 , which are the slowly varying components of 
mechanical motion x{t) around its resonance frequency SIm, 

x{t) ~ Xi (t) cos Qut + X2 (t) sin ^m* • (5) 

The error signal is obtained by mixing these quadratures at 
the mechanical resonance frequency Qm with tunable phase 
and amplitude (see Methods), and is fed back into the cavity 
via amplitude modulation the pump beam. 

Using this scheme we demonstrate accurate feedback con- 
trol of the fundamental out-of-plane mode by tuning $. Fig- 
ure 3(c) shows the feedback induced resonance frequency shift 
(Sloff — nM)/27r (left axis, blue dots) and normalized efi^ective 
damping Fcft/rM (right axis, red dots) versus the feedback 
phase $ with a constant g. The experimental data of Fig- 
ure 3(c) is fitted using Equations. |3] and |4] giving excellent 
agreement. In both cases we use the same fitting parameter 
g = 30.7, thus demonstrating the reliability of our scheme. 

We particularly investigate the case $ = it/2, which 
corresponds to purely dissipative feedback and has re- 
cently received increased interest as it enables to cool a 
mechanical degree of freedom [30l [41] |42l |44l |45] . We apply 
dissipative feedback with different gain amplitudes to the 
fundamental out-of-plane mode of the oscillator (cf. Fig. 
3(a)). In presence of feedback the oscillator is described 
by thermodynamical equilibrium as long as the spectral 
density of thermal motion S^xfO — f^ivi] is larger than the 
readout imprecision Si[fl ~ ^m]. In this limit an effective 
temperature T^s can be defined as T^g = (rM/rcff)r. Thus, 
starting from room temperature, we achieve a minimum 
effective temperature T^s — 700 mK, limited only by the 
readout sensitivity. To demonstrate the versatility of our 
method, we perform feedback cooling on the third order 
harmonic of the out-of-plane mode by simply changing 
the modulation/demodulation frequency. The results are 
shown in Fig. 3(d) and demonstrate that the scheme can be 
efficiently used for different modes. 

Force resolution enhancement via feedback control. 

High-Q nanomechanical oscillators, combined with high sen- 
sitivity readout, are ideally suited for measurements of 
weak forces |2l |9]. Their sensitivity is mainly limited by 
the thermal Langevin force Fth with the spectral density 
(SD) given by the fiuctuation-dissipation theorem S'ppffJ] = 
AksT /Qlm{l/x['^]), where fee is Boltzmann's constant. In 
the limit of high-QM harmonic oscillators, the SD reads 
Sppfn] ~ 4meftfeBT'rM and is proportional to both the mass 
and mechanical damping. For our integrated system this sen- 
sitivity is Sppfil] ~ (74aN)^/Hz, comparable to the best val- 
ues reported at room temperature for cantilever-based sys- 
tems [40, however with the benefit of an almost two orders 
of magnitude higher bandwidth [21]. In principle, detecting 
stationary incoherent forces embedded in the thermal noise is 



possible by means of energy averaging. However, this aver- 
aging converges very inefficiently, scaling with r^^* [T3| and 
thereby limiting the force detection threshold. In the follow- 
ing, we present a detailed study of energy averaging, and show 
that its convergence can be drastically improved by using dis- 
sipative feedback. 

We are interested in detecting a signal Xsig(t) driven by 
a stationary incoherent force 5Fsig{t), which is uncorrelated 
with the thermal force and for which S'ppfSl ~ Q^-^] <^ 
SppfrJ ~ Qm] holds true, i.e. the incoherent force is buried 
in the thermal noise. In presence of 5-Fsig, the energy of the 
transducer = ^mcs^mi^'^) can be written as the con- 
tribution of two independent terms Sth = ^nT-cB^i-nixti^) and 
iSsig = ^rncs^mi^sig) ) whcrc (...) denotes the average over the 
statistical domain describing a; as a random variable. Thus, 
comparing the transducer's energy with and without SF^ig be- 
ing applied enables in principle its detection. Experimentally, 
statistical averages are not available, and one has to rely on 
estimators. Assuming J-Fsig to be a Gaussian white noise for 
simplicity, s{t) = 1/r JJ^ dtX^{t) is a non-biased estimator for 
the measurement of the energy, with X denoting any of the 
motion's quadratures. X{t) is the sum of two uncorrelated 
contributions Xth(t) and Xsif,{t), giving therefore 

s(r) = Sth-th(r) + 2sth-sig(T) -|- Ssig-aig(r) (6) 

with Si_jg{th,sig} (■'") = I/t d,tXi{t)Xj{t). If we assume 
the measurement sampling time Ts to be small compared to 
the mechanical coherence time Ts ^ I/Fm, the statistical 
dispersion of each of these 3 terms can be calculated [47] to 
be: 

Asi_i(T) = V2H{rM,r){Xf), (7) 
Asth-sig(r) = H{rM,r){X^y,)i{x!,g)K (8) 

where //(FM,r) = ^2 (e^^"" - 1 + Fmt) ' / (Fmt) de- 
pends only on the mechanical linewidth of the transducer. 
Noting that {X^} ~ S'ppIQm] '^\xM\^ , it follows that 
{X^ig) <c: {X^^^), and we can write As(r) ~ Asth-th(T): Av- 
eraging X^ along r enables detecting a minimum excess of 
energy AS — |meffr2mAsth-th('r), related to a minimum de- 
tectable force with a SD Spp'™'"[f2m, t] given by: 

osig,minrQ i ASth — th('^) 

= V2H{rM,T)S%[ftM]. (9) 

Finally, we define the equivalent force spectral resolution 
(EFSR, in aN/\/Hz units) as the square root of the minimum 
detectable force SD : 

5F[Qm, r] = ^5;'|'"""[r!n„r]. (10) 

Equations [9] and [To] show that the time evolution of the EFSR 
is fully described by the function H{Tm,t). For r <^ I/Fm 
the EFSR stays constant due to the mechanical memory of the 
transducer: averaging over a time smaller than the mechani- 
cal coherence time is completely inefficient, corroborating the 
advantage of using high frequency, high-Qivi oscillators. For 
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T ^ I/Fm, the EFSR varies as (Fmt)~3 confirming previ- 
ous reports [13j . The discussion above shows that obtaining a 
given force resolution becomes faster as the damping rate of 
the oscillator increases. As dissipative feedback can provide 
an effective damping rate significantly larger than the intrin- 
sic one, while keeping the thermal force unchanged, it enables 
substantial decrease of the averaging time. In presence of dis- 
sipative feedback, the EFSR is 



\/rM/Feff X 5F[an,r], (11) 

where Fefi and SF[VLm,T\ are given by Equations [4] and [9] 
Consequently and importantly, feedback control reduces the 
time needed to reach a given force resolution by a factor of 
Feff/Fin which is of high practical relevance. 

Figure 4(a) shows the trajectory of mechanical motion 
in its phase-space, without and in presence of dissipative 
feedback (blue and purple traces, respectively). Figure 4(b) 
shows the EFSR calculated as a function of VbS and r using 
Equations [9] and 10 Figure 4(c) shows the corresponding ex- 
perimental results. The contours correspond to plane regions 
with constant resolutions described by (Foflr ~ ai)^.^^+. 
The good agreement between theory and experiment con- 
firms that an increase of Fch enables decreasing r by the 
same amount, while maintaining the same force resolution: 
detecting a weak stationary incoherent force is 57 times 
faster in presence of a feedback rate Fcft/F = 57 (upper 
dashed line in Figure 4(c)) than without feedback (lower 
dashed line in Figure 4(c)). Figure 4(d) shows details of 
the EFSR time evolution corresponding to these two cases. 
The upper curve corresponds to Fcfr/FM = 1 (no feedback) 
and the lower to Vcs/Vm = 57. Both curves show a scaling 
with r^^* in their asymptotic behavior (dashed gray lines), 
in excellent agreement with theoretical expectations. A 
curvature is present in the upper curve at short averaging 
times, which is absent in the lower one. We attribute this 
feature to the mechanical memory time of the transducer 
= I/Fm — 26 ms as described by the expression of 
//(FmjI") given above. Figure 4(e) shows the short term 
evolution of the force resolution in absence of feedback 
(framed region in Fig. 4(d)) together with our theoretical 



model (Eqs. [9] and 
parameter being usee 



T0[ see Methods) with no adjustable 
Excellent agreement is again found 
between measurement and theory, thus validating our model 
for weak force detection. Fig. 4(d) gives the improvement 
of the EFSR as a function of Feff/FM, for a given averaging 
time of r = 0.3 s as given along the vertical dashed line in 
Fig. 4(b). The straight line corresponds to a power-law fit 
with an exponent found to be 0.27, in agreement with a value 
of 1/4 expected from the theoretical discussion. 

Feedback assisted incoherent force detection. We now 

explicitly show how feedback allows to enhance the system's 
force resolution by detecting a weak, incoherent radiation 
pressure force. Particularly, we show that feedback assisted 
force detection may even allow the detection of incoherent 
signals that would remain out of reach under standard 
operating conditions. We apply the force to be detected by 
sending a third laser tone into the tapered fiber which has a 
randomly modulated amplitude (see Fig. 2(a) and Methods). 
It is important to emphasize that feedback enables an 
improved force resolution only for signals with coherence 



times being small compared to the effective mechanical decay 
time 27r/Fcfi (i.e. spectrally broad signals c.f. Fig. 5(a)). 
By measuring the autocorrelation function of our signal, we 
determined its coherence time to be more than three orders 
of magnitude smaller than the intrinsic mechanical decay 
time 27r/FM (see Fig. 5(b) and Methods), so that feedback 
gains up to a few hundreds can be safely employed. To 
detect the signal we use the "on/ofF' protocol suggested 
above: We first average the probe phase fiuctuations while 
the signal is being applied and compare these to the averaged 
probe fluctuations in absence of the signal. One critical 
limit of this method results from sensitivity fluctuations 
giving an upper boundary for the possible averaging time. 
With our system we determined the sensitivity drifts to 
be on the order of 1% at the minute scale: Signals with 
S^Kfi ~ Qm] < 0.01 X S%[Q. ~ ^m] can therefore only 
be detected within less than one minute of averaging using 
the on/off protocol. Reducing the required acquisition time 
via feedback is therefore crucial, as it enables breaking 
the sensitivity limits induced by the stability constraints 
of the system. We implemented the on/off protocol for 
various signal powers, feedback gains, and averaging times. 
Figure 5 shows signal detection for three different power 
levels, deduced from the signal's autocorrelation functions 
shown in Fig. 5(c). Figs. 5(d,e,f) show the estimated energy 
excess induced by the presence of the signal (s(t) — fth)/fth 
(expressed in units of thermal energy fth) as a function 
of averaging time r and feedback gain. The color of the 
points is related to the averaging time, varying between 
s and 35 s from blue to red. For each signal power the 
estimated energy excess converges towards the expected 
value (dashed blue lines, inferred from the level of the 
autocorrelation), with the final dispersion being in excellent 
agreement with our model (pink lines, see Methods). In 
particular, it is clearly confirmed that employing high 
feedback gains enables a significant decrease of the measure- 
ment imprecision for a given averaging time. Finally, we 
focus on the detection of the signal with the lowest power 



(5;'|[n ~ Hm] =^ 0.04 X Sf^[n ~ ^m] (15aN/VHz)^), with 
the results shown in Figs. 5(f) and 5(g). The pink region in 
Fig. 5(f) represents what we define as the Detection Zone 
(DZ) highlighting the points located within five standard 
deviations of the expected value of the measurement. The 
signal is considered to be detected in a given averaging time 
after which the centered normalized energy estimator (CNE) 
(s(r) — fth)/£th stays confined in the DZ. Fig.5(g) shows 
the detailed time evolution of the CNE given by the framed 
regions in Fig. 5(f). In absence of feedback (Fig.5(g)(i)), 35s 
of averaging is clearly insufficient to make a conclusion about 
the presence of the signal. Setting a moderate feedback gain 
g = 6.5 (Fig. 5(g) (ii)) reduces the measurement dispersion 
to 2 standard deviations after 35 s of averaging, but is still 
insufficient to decide upon signal detection. Fig.5(g)(iii) 
shows that employing a relatively high gain gr = 32 enables 
detecting the signal, i.e. the applied incoherent force, within 
the same 35 s of averaging duration. Note that reaching this 
resolution without the assistance of feedback would require 
averaging 32 times longer amounting to approximately 
20min. However, the sensitivity drifts in our system make 
averaging inefficient for this timescale. Thus, feedback is the 
condition sine qua non in this experiment. 

Conclusion. In conclusion, we have realized an integrated 
hybrid force transducer system for sensing applications of ul- 
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trasmall signals. The integration of the system is thereby of 
crucial importance, as it allows a stable operation over an ex- 
tended period of time. Otherwise, even the smallest changes 
in the system, particularly regarding the evanescent coupling 
between mechanical beam and the optical cavity, would pro- 
hibit the detection of small signals well below the thermal 
noise. In contrast to other high-Qivr and ultra-sensitive trans- 
ducers [5l[T6], our system operates at room temperature and 
profits from well-established read-out techniques at optical 
frequencies. Moreover, our system also strongly benefits from 
the optomechanical interaction, as it enables to use multiple 
optical modes as channels for interaction (readout and actua- 
tion) with the nanomechanical beam. Concurrently, the inte- 
grated system could be also implemented in a cryogenic envi- 
ronment in a straightforward manner |48| . The proposed and 
experimentally implemented scheme to improve the force sen- 
sitivity for incoherent signals via dissipative feedback has the 
potential to significantly reduce measurement times for detec- 
tion of weak signals. This scheme could be potentially applied 
for magnetic resonance force microscopy [13] and other signals 
such as transitions in nitrogen- vacancy (NV) centers [49ll50j . 
The system can be readily functionalised, for example via 
the use of nanomanipulators [51] or using a scanning atomic 
force microscope 52, 53 . Finally, the integrated system is 
a promising candidate for experiments in optomechanics ^36] 
due to its high Qm, potential of cryogenic implementation 
and applicability to feedback cooling methods that do not re- 
quire the resolved sideband regime for ground state cooling 
|44| . The lowest occupancy achieved with our system at room 
temperature is fi « 5000 (Teff = 700 mK) giving a promising 
outlook for operation at cryogenic temperatures. 

Methods 

Fabrication The samples are fabricated on a silicon wafer with 
thermally grown 2 fj,m silica and 100 nm LPCVD deposited high- 
stress stoichometric Si3N4 thin films on top. First, the nanobeam is 
created using electron-beam lithography followed by a timed C4F8- 
SFe chemistry reactive ion etching. In a second lithography step, 
the polymer mask for the disk resonator is defined in such a way 
that the disk partially overlaps the nanobeam. The disk resonator 
was created via a wet etch in a buffered hydrofluoric bath. Due 
to the isotropic nature of this etch, material below the polymer 
mask is partially etched resulting in an angled sidewall of the disk. 
Using this effect and ensuring the correct alignment of the disk 
mask with respect to the nanobeam, the release of the nanobeam 
is achieved with the structure being positioned above the wedge of 
the disk resonator. Release of the combined system is accomplished 
by anisotropic silicon etching in a potassium hydroxide bath. 

Actuation of the nanomechanical oscillator A pump-probe 
scheme is used to actuate and read out the motion of the nanome- 
chanical oscillator. Both the pump and the probe light fields are 
provided by external cavity diode lasers (ECDL). The probe beam 
is phase modulated using an electro-optical modulator (EOM) and 
the pump is amplitude modulated via a LiNbOs Mach-Zehnder 
modulator (AM). Both beams are combined with a fiber-based 
beam splitter (CI) and coupled to two different optical modes us- 
ing a tapered fiber. Both optical modes belong to the same WGM 
mode family and are separated by one free spectral range. The light 
coupled back from the resonator to the tapered fiber is split by a 
second fiber based beam splitter (C2) and tunable optical filters 
(TOFl and TOF2) are used to separate the spectral signatures of 
the pump and the probe beam. Light is detected via fast photo re- 



ceivers (PDl and PD2) and analyzed with an oscilloscope as well as 
an electrical spectrum analyzer (ESA) or network analyzer (NA). 
The frequency of both lasers is locked on the side of the fringe of 
the respective optical modes by providing an error signal to the 
piezo controller of the laser. In case of the ringdown measurement, 
the pump is amplitude modulated at the mechanical resonance fre- 
quency by providing to the AM the corresponding radio frequency 
signal from the NA set to zero span mode. The detected signal is 
analyzed by the NA. After interrupting the modulation a temporal 
decrease of the oscillation amplitude can be observed from which 
the mechanical Q is obtained. In case of the response measurement, 
the signal provided to the AM from the NA is swept around the 
mechanical resonance frequency and the response of the system 
is demodulated and recorded by the NA. The sweeps are always 
conducted from smaller to larger frequencies. 

Implementation of the feedback mechanism The experimen- 
tal setup is in principle the same pump-probe scheme as described 
above, however with the difference that the signal fed to the AM is 
derived from quadrature detection. Quadrature detection is accom- 
plished by splitting the signal from the PD detecting the light of the 
probe beam and demodulating it. The demodulation is achieved 
by mixing (low-noise analog mixer MINICIRCUITS ZP-3-S-I-) both 
arms of the splitted signal with a radio frequency source (at the 
mechanical resonance frequency Ojv//27r) from a signal generator 
(Tektronix AFG3102) with the phase of the source delivered to 
both arms being shifted by 90°. The corresponding signals are fil- 
tered with a low-pass (3 dB cutoff at 300 kHz) and give the time 
evolution of both quadratures. To obtain the error signal for the 
feedback, both low-pass filtered signals are modulated at the me- 
chanical resonance frequency by a second signal generator with the 
phase being again shifted by 90° for both arms. The phase between 
both signal generators is locked to a 10 MHz voltage controlled os- 
cillator. Adjusting the phase difference corresponds to adjusting 
the phase <& of the feedback signal. Both signal arms are combined 
after the modulation and after a low noise amplifier (Miteq 1447) 
fed to the AM. 

Enhancement of force resolution using feedback To de- 
termine the measurement accuracy (Eq. [tJ corresponding to 
a given feedback gain, we first extract the quadratures of mo- 
tion by sending the detected signal to a spectrum analyzer 
(Agilent MXA9020A), operating in I/Q mode with a sampling 
rate of I/ts ~ 7 kHz. We record independent realizations of 
the thermal noise {Xi{t), X2{t))ij along A'^ different acquisition 
times Tj ranging from 27 ms to 3.77 s. Each of the resulting se- 
quences is used in order to determine the corresponding averaged 
energy s,j(t,) = I/tj f,J' dt (^^ (^^^.,(4) + (^,^.) (t)) • The mea- 
surement resolution after averaging time tj is then determined as 
As(tj) = ^Var ■ We define the equivalent force resolu- 

tion (represented in Fig. 4) as the equivalent motion resolution nor- 
malized to the area of the mechanical amplitude response 5F{Tj) = 

J As{Tj)/ If |x[f^]P> the latter being determined by measur- 
ing the expected value of the energy (s(Tj)) ~ 1/A'" y^j^p s(t^), 
using that {s{r)) = [Hm] J^T '^IxMl^- 

Feedback assisted force detection To obtain the incoherent 
radiation pressure force, the output laser field of a third ECDL 
is sent into an AM fed with a Gaussian white noise generated 
by a function generator (Agilent 33250A). To avoid saturating 
the AM, the generator output is filtered using a selective 3 MHz 
bandwidth low-pass filter. The optical frequency of this ECDL is 
tuned in order to match a third optical resonance of the WGM 



6 



cavity, distinct from the one used to probe (ECDLl) and to ac- 
tuate (ECDL2) the feedback response. This third "noisy" tone 
is sent into the tapered fiber using a fiber based beam spHtter 
(C3 on Fig. 2(a)) and is later separated from the other two tones 
using the combination of another beam splitter (C4) and of an 
optical filter (TOF3). A photodetector (PD3) enables measur- 
ing the transmitted intensity fluctuations. To determine the au- 
tocorrelations functions presented in Fig. 5(b), we proceeded in 
the following manner: The time evolution of the signal's quadra- 
tures is first recorded for a duration Ts ~ 35 s at a sampling rate 
dt ~ 0.2 ms using an ESA operating in I/Q mode. The autocor- 
relation being quadrature-invariant for Gaussian signals, only one 
quadrature X is conserved for its computation, which we restricted 
to delay times St smaller than (5Tmax = 30 ms. For any 5t £ 
[— <5Tmax, <5Tmax] fl dt X N the autocorrelation value {X{0)X(St)) is 
obtained by evaluating the discrete time integral {X(0)X(St)) ~ 

J, /ir^"''^""'^ AtX{i')X(t-\-&T'), the latter equality applying 

by virtue of the ergodic theorem. As both signals being measured 
(the intensity noise and the thermal displacement of the nanome- 
chanical oscillator) are resulting from white noises, their autocor- 



relation function is expected to be the impulse response of the 
measurement apparatus used to detect them ({nanomechanical os- 
cillator-l-WGM resonator-|-photodctector-|-ESA} and {WGM res- 
onator+photodetector-|-ESA} for the mechanical motion and the 
intensity noise, respectively). For the probe signal, this impulse 
response is mostly due to the mechanical response and can be ap- 
proximated by a decaying exponential, since the sampling rate is 
dt <C 27r/rM. For the intensity noise, this impulse response promi- 
nently involves the 5 kHz wide frequency gate set by the I/Q mode 
parameters of the ESA, such that the autocorrelation linewidth in 
Figs. 5(b,c) is not related to the coherence time of the intensity 
noise but simply to this frequency gate. The pink dispersion lines in 
Fig. 5(f) are determined using that A£(t) = Imefff^m^'^th— th(^) 
since the low SNR hypothesis applies. This is, however, not the 
case in Figs. 5(d,e), where the SNR are (Argjg)/(X^j^) ~ 5 and 
(Jf^jg)/(X^j^) ~ 0.5 respectively. The corresponding dispersions arc 
therefore evaluated using an extended model taking into account 
the terms Assig_tii(''') and Assig_sig(T) (Eqs. [t] and [sjl . 
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Figure 1: Hybrid nanomechanical transducer system on a chip, (a) Scanning electron image (false colors) of the hybrid 
on-chip system consisting of a doubly-clamped, high-stress Si3N4 nanomechanical beam (blue) near-field coupled to a silica disk 
resonator (red). Both, the resonator as well as the pads holding the beam rest on silicon pedestals (green), (b) Magnified image 
of the proximity area, (c) Optical micrograph of a tapered fiber coupled to the hybrid system, (d) Finite-element simulation 
of the electromagnetic field distribution of the fundamental optical mode confined in the disk resonator. The nanomechanical 
beam is coupled to the evanescent near-field of the optical mode, (e) Room-temperature Brownian noise of a nanomechanical 
beam with a fundamental resonance frequency of 2.88 MHz and rrics = 9 pg. Fitting the mode gives a mechanical quality 
factor of Qm ~ 4.3 x 10^. Inset: Finite-element simulation of the beam's fundamental out-of-plane mode. 
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Figure 2: Actuation of the nanomechanical oscillator by radiation pressure, (a) Scheme of the actuation-readout 
setup used for the actuation and feedback experiments. By amplitude modulating a third laser (green shaded area) an incoherent 
force can be applied on the transducer for direct force detection experiments. Details of the setup are given in Methods, (b) 
Ringdown measurement of the fundamental mode fitted with an exponential decay function (red line) confirming the high Qm • 
(c) Response measurements for an increasing actuation of the nanomechanical motion. The actuation power is increased by 5 
dB for the curves from bottom to top. The direction of the sweep is from lower to higher frequencies. With increasing actuation 
the response of the oscillator is increased and above a certain threshold a mechanical bistability sets in. The middle curve 
(dark yellow) was taken at the threshold and is used to calibrate the actuation in terms of the critical displacement Ac. Oc is 
found to be 1.6 nm for our system. 
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Figure 3: Feedback control of the nanomechanical transducer, (a) Feedback cooling of tiie fundamental out-of-plane 
mode. The light blue curve corresponds to a squashing of the background present. While it can still be fitted our model, we 
do not use it to determine Tes. Inset: Finite-element simulation of the beam's fundamental out-of-plane mode, (b) Feedback 
scheme based on demodulation and modulation of the incoming probe signal (Details see text and Methods), (c) Feedback 
cooling of the fundamental in-plane mode. Inset: Finite-element simulation of the beam's fundamental in-plane mode, (d) 
Feedback cooling of the third order harmonic of the out-of-plane mode family. Inset: Finite-element simulation of the beam's 
third order out-of-plane mode, (e) Frequency detuning (Sleff — ^^m) /27r of the fundamental out-of-plane mode (left axis) and 
ratio of the effective linewidth to the intrinsic linewidth (right axis) versus the phase difference $ between the demodulation 
and modulation signals. 
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Figure 4: Force resolution enhancement via feedback control, (a) Evolution of mechanical motion in phase-space 
without and in presence of feedback (blue trace and purple traces, respectively). The histograms correspond to the statistical 
distributions of quadratures Xi and X2. (b) Calculated force resolution as a function of the effective damping TcS and average 
time T and (c) corresponding experimental result, (d) Evolution of the force resolution as a function of r without feedback 
(upper curve, corresponding to the blue dashed line in Fig. 5(c)) and with Feff/FM = 57 (lower curve, corresponding to 
the purple dashed line in Fig. 5(c)). Both curves show very good agreement with an asymptotic behavior scaling with r^'''' 
(dashed gray lines), (e) Short term evolution corresponding to the framed region in (d) together with the theoretical model, 
(f) Evolution of force resolution as a function of Foff/FM for a fixed averaging time of t = 0.3 s, corresponding to the evolution 
along the vertical dashed line in Fig. 5(b). 
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Figure 5: Feedback-assisted force detection, (a) Scheme illustrating the validity limit for the feedback-assisted weak 
random force detection. Left: The quadrature spectrum of the driving signal (in blue) has to be broad when compared to the 
mechanical spectrum (in purple). Right: Equivalently, the autocorrelation function of the signal (linked to its spectrum via a 
Fourier Transform (FT) operation) has to be narrow as compared to the autocorrelation function of the mechanical motion, (b) 
Experimentally measured autocorrelation functions associated with the signal intensity noise (in blue) and with the thermally 
driven mechanical motion (in purple), (c) From top to bottom: Intensity noise autocorrelation function for decreasing noise 
powers, associated with the signals detected in (d), (e) and (f) respectively, (d-f) Signal contribution to the mechanical energy 
(expressed in units of thermal noise), measured along various averaging times and feedback gains. Each point corresponds 
to the signal contribution obtained in presence of a feedback gain given by its abscise value and after a given averaging time 
related to its color (single shot measurement for the bluer poins, 35s of averaging for the redder ones). The blue dashed lines 
correspond to the expected values of the signal contributions to the mechanical energy (5.25 x £th, 0.5 x fth and 0.04 x fth, 
respectively), deduced from the levels of the signal autocorrelation functions shown in Fig. 6(c). The pink lines correspond 
to the gain evolution of the energy dispersion expected after 35 s of averaging (see Methods). The pink region in (f) gives the 
signal detection zone, which we define as corresponding to a measurement accuracy of 5 standard deviations or better, (g) Figs. 
6(g) (i), 6(g) (ii) and 6(g) (iii) represent the detailed evolution of the signal contribution with averaging time for feedback gains 
of 0, 6.5 and 32 respectively for the lowest signal applied, E^ig — 0.04 fth. The blue dashed line corresponds to the expected 
measurement dispersion (see Methods). The pink region corresponds to the detection zone. 



